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A Maximum Stellar Surface Density in Dense Stellar Systems 
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ABSTRACT 

We compile observations of the surface mass density profiles of dense stellar systems, 
including globular clusters in the Milky Way and nearby galaxies, massive star clusters in 
nearby starbursts, nuclear star clusters in dwarf spheroidals and late-type disks, ultra-compact 
dwarfs, and galaxy spheroids spanning the range from low-mass, "cusp" bulges and ellipticals 
to massive "core" ellipticals. We show that in all cases the maximum stellar surface density 
attained in the central regions of these systems is similar, S^ax ^ 10"MQkpc^^ 20 g 
cm~^), despite the fact that the systems span ^ 7 orders of magnitude in total stellar mass 
M*, ^ 5 in effective radius R,,, and have a wide range in effective surface density M^/R^. The 
surface density limit is reached on a wide variety of physical scales in different systems and is 
thus not a limit on three-dimensional stellar density. Given the very different formation mech- 
anisms involved in these different classes of objects, we argue that a single piece of physics 
likely determines E^ax- The radiation fields and winds produced by massive stars can have a 
significant influence on the formation of both star clusters and galaxies, while neither super- 
novae nor black hole accretion are important in star cluster formation. We thus conclude that 
feedback from massive stars likely accounts for the observed Emax, plausibly because star for- 
mation reaches an Eddington-like flux that regulates the growth of these diverse systems. This 
suggests that current models of galaxy formation, which focus on feedback from supernovae 
and active galactic nuclei, are missing a crucial ingredient. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — star formation: 
general — cosmology: theory 



1 INTRODUCTION 

Feedback from massive stars and black holes is widely believed to 
play a critical role in the formation of other stars, stellar clusters, 
and entire galaxies. However, the precise physical mechanism(s) 
that dominate the coupling between stars and/or black holes and 
their surrounding environment have not been definitively estab- 
lished. Unlike the energy deposited by, e.g., a jet, the momen- 
tum produced by a central star or accreting black hole cannot be 
radiated away. When the momentum deposition rate exceeds the 
strength of gravity, gas can be efficiently expelled from the system. 
This limit is analogous to the Eddington limit from stellar physics; 
however, the gas in star-forming regions and galaxies is dusty, and 
thus the appropriate opacity is not that of electron scattering, but 
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the much larger opacity due to dust (see e.g. IScoville et alJIlOOll : 
iMurrav et alJlOOalThompson et al]|2005l) . 

This generalization of the Eddington limit has been posited as 
an explanation for the characteristic sizes of massive stellar clusters 
( lMurravll200^ , as well as the maximum luminosities of the most 
rapidly star-forming objects in the Universe, bright sub-millimeter 
galaxies with ~ lO'^ Lq inside of ^kpc (star formation rates > 
1000 MQyr"' kpc"^ e.g.. Younger et al.l2008l:IWalter et alj2009h . 
Outflows driven by radiation pressure may help explain phenom- 
ena as diverse as turbulence in star-forming regions, the relation 
between galaxy mass and stellar population metallicity, the corre- 
lations between the masses of black holes and their h ost spheroids, 
and the abundance patterns of the IGM (see e.g. iMurrav et al.l 
2005"; 'Thompson et alj |2005l: lOppenheimer & P avel l200 d l2008l: 
Ki-umholz et al. 2009; IMurrav et alj l2009l : iKr u mholz & Matznei) 
2009). Small changes in the nature of stellar feedback — for ex- 
ample, whether or not winds are really momentum or pressure 
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driven, or the characteristic scales of wind driving — have dra- 
matic consequences for predictions of models of galaxy forma- 
tion, affecting the global stru cture of even Mi lky Way-mass sys- 
tems JOovemato et alJlOOTl i lScannapieco et alJboOSh . 

Out flows appear to b e ubiquitous from rapid ly star-forming 
galax i es iHeckinan et aiT iggO: Martin 1999, 2005; Heckman etal] 
120001 ; lErb et"ai]|200d IWeiner et al. 2009). However, quantifying 
the basic properties of these outflows is challenging; their mass- 
loading, dynamical importance, and implications for future star for- 
mation remain subjects of considerable debate. This uncertainty 
translates into uncertainty in the dominant physical processes driv- 
ing such outflows (e.g., radiation pressure, supernovae, AGN vs. 
star formation, etc.). It is also unclear whether or not feedback pro- 
cesses and the outflows they generate actually leave an imprint on 
the structural properties of star clusters and galaxies. 

In this Letter, we show that a diverse set of dense stellar sys- 
tems — from star clusters to the centers of elliptical galaxies — 
have a maximum stellar surface density. We argue that this obser- 
vational fact provides a strong constraint on models of feedback 
in star cluster and galaxy formation, pointing to the critical role of 
massive stars, rather than supernovae or AGN. In § [2l we present 
our compilation of observations of dense stellar systems (Table[Tll, 
and compare their surface density profiles and maximum stellar 
surface densities. In § [3] we discuss the implications of our re- 
sults and interpret the observed maximum stellar surface density 
in terms of models in which star formation reaches the (dust) Ed- 
dington limit. 

Through out, we adopt a JIm = 0.3, SIa ~ 0.7, h = 0.7 cos- 
mology and a IChabriejJiooah stellar IMF, but these choices do not 
affect our conclusions. Changes in the IMF would systematically 
shift the stellar masses of the massive systems, but reasonable (fac- 
tor < 2) changes are within the scatter in the data we compile. 

2 OBSERVATIONAL RESULTS 

In Table [T] we summarize our compilation of a large sample of 
the observed stellar mass surface density profiles of dense stellar 
systems, including globular clusters (GCs) in the Milky Way and 
nearby galaxiesQ the nuclear stellar disk observed around Sgr A* 
in the galactic center, massive star clusters in starburst regions, nu- 
clear star clusters in dwarf spheroidals and late-type disks, ultra- 
compact dwarfs, and classical spheroids spanning the range from 
low-mass, "cusp" bulges and ellipticals (those with steep central 
surface brightness slopes) to massive "core" galaxies (these pre- 
dominate at the highest masses, with shallower central slopes). 
We also include massive compact ellip ticals observed at h igh red- 
shifts, z ~ 2 — 3; as a lready shown in iHopkins et al. I ( l2009a) and 
ISezanson et al.l j2009l) . these compact high z systems in fact have 
the same maximum density as the central regions of today's ellipti- 
cals. 

For each object, we adopt the best-fit surfa ce brightness pro- 
file determined by the authors. This is us ually alKing JT966I) pro- 
file for the less massive systems, and a ISersid ri968l) or Nuker 
profile for the more massive objects. For the MW nuclear stellar 
disk we adopt the total mass from Lu et al. ( 2009) and the pro- 
file shape fitted in IPaumard eTal] ilOOd) (E, cx (/?^ + i^Le)"')l3 



' We do not distinguish between GCs that might have undergone core col- 
lapse and those that have not (see §[3). 

^ We use the notation S«, Sgas, and Stot to refer to the stellar, gas, and total 
(stellar plus gas plus dark matter) surface densities. Throughout, the terms 
gas mass or gas surface density refer specifically to "cold" (i.e. star-forming 
or rotationally supported disk) gas. 



Table 1 . Observations of Dense Stellar Systems 



Reference ' 


Description" 


Symbol^ 


Lu et al. (2009) 


MW Nuclear Disk 


dark green square 


^arris(1996) 


MW GCs 


pink circle 


Barmby et al. (2007) 


M3 1 GCs 


yellow circle 


^eikuba et al. (2007) 


Massive Cen A GCs 


green pentagon 


McCradv & Graham ^2007) 


M82 SSCs 


orange X 


^akher et al. (2005) 


Sd+dSph nucleai- SCs 


red + 


Boker etal. (2004) 


Sd disk nuclear SCs 


magenta asterisk 


Geha et al. (2002) 


dSph nuclei 


violent inv. triangle 


^a§egan et al. (2005) 


M87 UCDs 


light blue star 


Evsti.tineeva et aL ^2007) 


Virgo UCDs 


dark blue triangle 


^ilker et al. (2007) 


Fornax UCDs 


cyan diamond 


jCormendv et al. (2009) 


Virgo Es 


red circle 


Lauer et al. (2007) 


local massive Es 


violet square 


van Dokkum et al. (2008) 


z — 2 — 3 massive Es 


orange star 



For most of the objects we adopt the fitted functional forms for the density profiles 
presented in each paper. For thelK ormend y et all 12009!) and'Lauer et al,' '2007^ sam- 
ples, however, the data are sufficiently high dynamic range that we use the observed 
profiles directly. Note that the profiles shown in . Lauer et ah i 2007) include only the 
HST data; the profiles used here are the composite (HST+ground-based) profiles used 
therein to determine effective radii. 

^GC^Globular cluster, SC^Star cl uster, SSC^ Super star cluster, UCD^Ultra- 
compact dwaif, E— elliptical. Note that lGeha et alj |20Q^ and jWalcher et alj |2005l l 
refer to their dwarf galaxy hosts of nuclear star clusters as dEs; we follow 
iKormendy et alj|2009| ) and classify these as dwarf spheroidals (dSph). 

^Plotting symbol used for objects in each sample throughout. 



The data from lKonnendv etaL I ( I2OO9I) and lLauer et all ( |2007|) span 
sufficient dynamic range that we directly adopt the observed (PSF 
de-convolved) profiles rather than a fit to the data; for the other 
systems in Table [T] the data do not span sufficient dynamic range 
to see significant deviations from the fits used here. In all cases, 
we convert the measured light profiles to stellar mass density pro- 
files by first adopting the authors' determination of total stellar 
masses (re-normalized to our Chabrier 2003 IMF where appropri- 
ate). We then assume a radius-independent stellar mass-to-light ra- 
tio to convert the light profiles to mass profiles. The assumption 
of a radius-independent M, /L is reasonable for the well-resolved 
galactic systems (nuclear SCs and Es), which are relatively old 
and have weak obse r ved color or stellar population gradients (see 
e.g.lCote etalj|2006l ; ISanchez-Blazguez et alll2007l ; iHopkins et d] 
|2()08[) . The lower-mass systems (e.g. GCs) are approximately single 
stellar populations, so this should be a reasonable first approxima- 
tion in those cases as well. In the high z ellipticals the central den- 
sities are unresolved; we thus extrapolate the best-fit Sersic profiles 
inwards. For low-redshift ellipticals, this extrapolatio n is typically a 
factor ~ 2 higher than the true central densities (see IHopkins et al.] 
l2009al for more details). 

The left panel of Figure [T] shows the effective surface stellar 
mass densities Mi^/nR' of these dense stellar systems as a function 
of their total stellar mass, where Re is the projected half-stellar mass 
radius. The wide range of systems shown in Figure [T] differ sig- 
nificantly in their large-scale structure and formation mechanisms. 
Various works have compared these different object classes in this 
and other projections of the fundamental plane, e.g., in the space 
of effective surface or luminos ity density, mass, mass-to-li^ht ra- 
tio, size, or velocity dispersion (Kormendy'l985';'Geha et al."200i 
Hasegan et al. 2005; Rejkuba et al. 2007; Evstigneeva et al. 2007|; 
[Pabringhausen et alj2008l; IHopkins et al.i2009bD . In some of these 
correlations there is apparent continuity from star clusters to mas- 
sive galaxies, but in others, different objects trace nearly perpendic- 
ular correlations and thus appear physically quite distinct. In terms 
of their effective densities, Figure[T]shows that massive galaxies are 
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Figure 1. Left: Effective surface stellar mass density Seff = M* (< Re)/iTR^ as a function of total stellar mass for the different samples described in the text 
and legend. Right: Maximum surface stellar mass density Smax obtained at small radii in the surface density profile (see Fig.[2j, as a function of total stellar 
mass, for the different samples described in the text and legend. The effective densities vary significantly across the different samples. However, there appears 
to be a well-defined upper Hmit ~ lO" MQkpc^^ to the maximum suiface stellar mass density; this limit holds for the full range of dense stellar systems 
considered here. 



significantly less dense than low-mass stellar clusters, and there is 
no indication of a universal maximum stellar surface density. 

However, E, is of course higher at small R<^Re. Figure[2]thus 
directly compares the surface stellar mass density profiles of these 
different systems. For clarity, rather than plotting every individual 
profile, we plot the median and ±1(t scatter in at each R 

for each class of objects listed in Table[T](we neglect the low-mass 
GCs, which have much lower densities in Figure [T] and remain rel- 
atively low-d ensity at all radii). Since the elliptical samples, in par- 
ticular that of lKormendv et alj ( l2009l) . span a very large range in Re 
and M*, we split them into 3 classes; low-mass ellipticals in Virgo 
(M, < IO'^Mq) and high mass cusp and core ellipticals. 

Figure |2] shows that the stellar mass surface density profiles 
appear to asymptote to a maximum Emax ~ 10"MQk;pc~^ ~ 20 
g cm^^, independent of object class, despite spanning an enor- 
mous range of M, and Re- The most massive Es maintain this 
density out to few 100 pc scales, with E, decreasing only weakly 
with radius. Massive GCs, on the other hand, also reach this sur- 
face density at small radii, but their densities then fall off rapidly 
with radius. Note that radii <Cpc are not resolved and are thus 
not plotted for these systems; extrapolations of the King model 
fits to arbitrarily small R, however, do not significantly exceed 
the maximum densities shown in Figure |2] Even the MW nu- 
clear disk, with its proxim ity to the nuclear black hole and evi- 
dence for a top-heavy IMF l lNavakshin & Sunvaevll200^ . does not 
appear to exceed Smax- Interestingly, recent observations suggest 
that there is a nuclear star cluster in the center of the Milky Way. 
The observations favor a radius-independent stellar surfac e density 
» 1.5 -3 X IO"MQkpc"- from ~0.02- I.Opc tSchoedel et all 
I2OO9I) ; this is comparable to E^ax found here. 

The right panel of Figure [T] shows the maximum stellar sur- 
face density E,„ax of each object as a function of its total stellar 
mass. For the King profile objects, Emax is defined to be the stel- 
lar surface density inside the fitted "core" radius (~ Rf). For the 
ellipticals, we define Emax as the (E,) inside the radius where E* 



falls a factor ~ 2 below the maximum E* measured[3 Reasonable 
variations in our definition of Emax make little difference to our 
general conclusion; although the effective densities of dense stellar 
systems (left panel) vary widely, many such systems (and most of 
the object classes) have roughly the same maximum stellar surface 
density ~ IO"MQkpc^^ at small radii, over a factor of ~ lO' in 
M*. 

3 DISCUSSION AND INTERPRETATION 

The fact that a wide range of star clusters, dwarf galaxies, and 
massive galaxies have the same maximum stellar surface density 
Emax ~ 10" Mq kpc~^ ~ 20 g cm^^ suggests that a common phys- 
ical process operated during the formation of these diverse systems 
to limit the stellar density that could be attained. It is worth em- 
phasizing that the observed maximum is a maximum surface den- 
sity. The different radii and masses sampled here imply that the 
same stellar surface density corresponds to very different three- 
dimensional stellar mass densities p. Whatever process regulates 
the growth of these systems selects a particular surface density, not 
a particular three-dimensional density. 

It is, in principle, possible that the observed Emax is due to the 
absence of sufficient gas reaching surface densities 2> Emax- For 
example, in galaxy mergers, the high star formation efficiency in 
pre-merger disks with high gas fractions and masses, as well as the 
inefficient angular momentum loss during gas-rich mergers, mean 
that it is difficult to generate gas surface de nsities far in excess of 
Emax on scales > 100 pc (see e.g. lHopkins et al. 2009c); this same 
conclusion may not, however, hold on smaller scales where E, is 

^ For core ellipticals, the maximum is well-defined. We caution, however, 
that for cusp ellipticals, it is not yet established whether or not the profile 
continues to rise at sub-pc radii. At these radii, though, the profiles are not 
necessarily representative of those at formation, as W-body relaxation and 
core-collapse can occur. We therefore focus on the maxima at the range of 
fixed small radii observed (~ 1 — lOpc), which are consistent (see Figure|2). 
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□ MW Nuclear Disk 




R [kpc] 

Figure 2. Stellar mass surface density profiles for the different samples. For each object class, thick solid lines (with points) show the median stellar mass 
surface density profile; shaded range shows the ±1 ct dispersion. The upper limit ~ lO" Mq kpc~^ is apparent, across a variety of classes. 



nonetheless observed to be < Smax (Fig.O. In the case of star clus- 
ters, the masses of the most massive giant molecular cloud (GMC) 
complexes present a similar limitation; on ~ 1 — lOpc scales, gas 
surface densities ^ E^ax would require GMC masses in excess of 
~ 10^ — 10^ Mq, at the limit of those inferred observationally. 

The observed maximum surface density occurs, however, over 
a huge dynamic range in mass (~ 10* — IO'^Mq), effective radius 
(~ 3pc — 50kpc), stellar population age (young clusters with ages 
^ lO^yr to spheroids and globular clusters with ages > lOGyr) 
and other properties. The formation mechanisms and formation 
timescales of these different classes of objects are also very dif- 
ferent. Whatever determines their maximum stellar surface density 
cannot, therefore, be specific to details of their formation or their 
global properties (e.g., M,, Re, etc.). As a result, it requires signif- 
icant fine-tuning for the same stellar Emax to apply independently, 
due to the availability of gas, in such a wide range of systems. In- 
stead, it appears likely that there is some rather generic physics 
that sets Emax- The key traits the systems we consider share are 
that they are baryon-dominated and likely formed in dissipational 
(rapid, gas-dominated) events. 

Because star clusters form and disrupt on a timescale much 
less than the time f or massive stars to explode as supernovae (e.g., 
lMurravetal]|2009l) , and because the time for stars to form at a 
gas surface density Egas ~ Emax is also less than the time for mas- 
sive stars to explode|j it is unlikely that supernova feedback is im- 
portant in setting Emax- In addition, the energy input from super- 
novae scales as £sn ~ (lO'^ergMQ ) E,i?", while the binding en- 

Note that the dynamical timescale is t^y^ ~ (GS^ax/'')^''^" ~ 1-5 X 
lO^rj^^ yr, where rio = r/ lOpc. 



ergy of gas scales as Ebini ~ GM^ /R ^ GEgasEtot/^^^- Although en- 
ergy feedback can in principle pick out a specific Emax, it is only 
below such a surface density, not above it, that feedback from en- 
ergy deposition is important; this is contrary to the phenomenology 
suggested by the data. 

Similarly, there is no evidence for energetically important 
black hole accretion during star cluster formation. The similarity 
of Emax in both star clusters and galaxies thus suggests that mas- 
sive stars themselves (on or near the main sequence) are critical for 
setting the observed Emax in dense stellar systems. Although we do 
not have a fully satisfactory theoretical explanation for Emax in the 
context of this interpretation, we briefly describe the possibility that 
it is set when these systems reach the dust Eddington limit. 

The Eddington luminosity (Lem) or flux (FEdd) for dusty gas is 

^Edd ftdd ^TtGc -1-1 ,,, 

— = — = PS 2500 K,o ergs s g , (1) 

M,ot Etot K.F 

where «io is the opacity in units of 10 cm" g^ ' , Mm and Etot are the 
total enclosed total (gas, stars, and dark matter) mass and surface 
density, and Kf is the flux-mean dust opacity. When the medium 
is optically-thick to the re-radiated FIR emission (Egas > 0.5 cm" 
g~') Kf can be approximated by the Rosseland-mean dust opacity, 
which has the form ac« ~ kqT' for T < 200 K and ~ 5 — 10 cm^ 
g~' ^ constant for 200 K <T< T^ub, where kq ~ 2 x 10^"* cm^ 
and Tsub ~ 1 500 K is the sublimation temperature of dust 
Semenov et al ][200i. The surface densities in Figures [T]&|2] are 
sufficiently high that if the systems had gas fractions of order fg ~ 
0.01 or larger during formation, then the medium was optically- 
thick in the FIR («K/gEmax > 1). For gas fractions of order unity, 
the temperature deep inside the optically thick gas was T > 200 K 
and kr ~ 5 — 10 cm^ g^' (for solar metallicity and Galactic gas-to- 
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dust ratio). These considerations motivate our scaling for n in the 
last equality of equation ([TJ. 

We can compare the Eddington flux to the flux from star for- 
mation implied by a given gas surface density if we extrapolat e 
the observed Kennicutt-Schmidt relation. From lKennicut 
E* ^ 2.5 X lO"''M0yr"' kpc"' [E.as /Mq pc"^]' '; other studies 
of high surface density systems suggest the index may be somewhat 
steeper, ~ 1.7 (Bouche et al. 2007). Together with the empirically 
calibrated relation between SFR and bolometric or total infrared lu- 
minosity, LsF = 2.22 X 1O^*W(M,/M0yr"'), this yields a simple 
scaling of flux from young stars versus surface gas density; 

— — « 2.8 ergs s g — (2) 

Egas VlO«MQkpc-2/ 

This relation shows that the flux from star formation reaches and 
then exceeds the threshold set by equation ([TJ for Egas ^ lO" — 
lO'^MQkpc"^ (depending on the index adopted), suggestively 
similar to the maximum we find. The Eddington flux can indeed 
be produced by a young stellar population (< lO'' * yr), for which 
the stellar light to mass ratio is v|/ = F* /S, ~ 3000 ergs s~ g^ 
for a standard IMF. A young stellar population is thus just capable 
of reaching the dust Eddington limit and supporting an optically- 
thick self-gravitating medium with radiation pressure. In this limit, 
however, it is unclear whether stellar radiation pressure can account 
for a maximum surface stellar density E^ax since both F« cx E* and 

Fedd OC Etot. 

Another possibility is that the relevant luminosity could, for a 
short period of time, come directly from the inflow of gas required 
to form the system, rather than from starlight. If a gas-dominated 
system of mass M and gas surface density Egas ~ Etot collapses on 
a free-fall time, the luminosity due to the release of gravitational 
binding energy is L ~ Mvjf oc Ega^^M''''', where M ~ M/[R/vff] 
is the inflow rate and v// is the free-fall velocity. Comparing the 
inflow limiinosity with the Eddington luminosity, it is straightfor- 
ward to show that there is a critical gas surface density above 
which the inflow luminosity exceeds the Eddington luminosity: 
E„it ~ 3 X lO"M0kpc~^/t;"o''^^M8^^ where Mg ^M/IO'^Mq is 
the mass at Egas ~ Ecrit, which is not necessarily equal to the to- 
tal mass of the system. This estimate is intriguing both because it 
is close to our observationally-inferred maximum surface density 
and because it provides a clear mechanism for obtaining a max- 
imum in the gas (and ultimately stellar) surface density: the flux 
produced by inflow oc Egas^ increases more rapidly with gas surface 
density than the Eddington flux oc Etot so the system would likely 
adjust to have Egas ~ Ecrit. However, the required inflow rates are 
extremely large, M ~ 2OOOM8''''M0 yr"' at Egas ~ lO" Mokpc"", 
and must be present even at very small radii (~ pc for star clusters 
and < 10 — 100 pc for massive galaxies). These required inflow 
luminosities are significantly larger than we have found in a pre- 
liminary analysis of our numerical simulations of the mergers of 
very massive gas-rich galaxies. As a result, it appears unlikely that 
the high inflow luminosities required to reach the dust Eddington 
luminosity can in fact be realized in all of the systems considered 
in Figure|2] 

One concern about invoking the dust Eddington limit is that 
the GCs in Figures [T] & [2] have low metallicities, typically ~ 0.1 
solar, but sometimes even lower. As a result the dust opacity will 
be ~ 10 times smaller than that used in equation (TJ and the massive 
stars present during GC formation would not reach ~ Fmd- We note, 
however, that the winds from massive st ars have a momentum flux 
comparable to that of the photons (e.g.. lLeitherer et al]|l999h and 
could have an impact similar to that attributed to photons here. 



Another concern about our conclusions drawn here is that 
the profiles observed now may not reflect those at formation, be- 
cause of A'-body relaxation. The relaxation time for a virialized sys- 
tem near E„,ax is Ir ~ lO'"^ yr [Etot/Emax]'''-[i?/10pc]'''l Note the 
steep dependence on R; for GCs, it is well-known that tn < lO'yr, 
implying they may be significantly affected by relaxation; but re- 
laxation cannot explain the > lOpc scales on which the centers 
of galaxies, spheroidal nuclei, and UCDs reach E„,ax. Caution is 
needed, however, at sub-pc scales, where the relaxation time is very 
short (~ 10** yr). Even in such cases, however, it is unlikely that the 
somewhat larger scales co nsidered here are strongly modified. Sim- 
ulations of core coflapse ( Cohn " 1 980l ; Ichemoff & WeinberdI 1 990l) 
tend toward sharp nuclear (<Cpc) mass concentrations and stellar 
mass segregation, but the mass fraction shifted is small, and the 
surface mass density on scales ~pc (the range of interest here) is 
largely unaffected (changing by a factor of < 2 — 3, within the scat- 
ter seen here). 

The maximum stellar surface density found here could readily 
be exceeded if gas "trickled in" slowly from large radii, and was 
allowed to form stars for a long time and at a low rate. In principle, 
there should be no limit to how large a stellar surface density might 
be attained in this case because all forms of feedback would be inef- 
fective at low star formation rates; in particular, the system would 
remain well below the Eddington limit. The interpretation given 
here requires that the high-density portions of the galaxy or cluster 
form in no more than a few, rapid events; this appears consistent 
with other observational indicators (e.g. the stellar populations). 
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